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V. D. Ivashchuk and V. N. Melnikov 

ABSTRACT 

A cosmological model describing the evolution of n Einstein spaces (n > 1) with m- 
component perfect-fluid matter is considered. When all spaces are Ricci-flat and for any a-th 
component the pressures in all spaces are proportional to the density: pf' = (1 — h^)p^ a \ 
= const, the Einstein and Wheeler-DeWitt equations are integrated in the cases: i) 
m — 1, for all h[ , ii) m > 1, for some special sets of h\ . For m = 1 the quantum worm- 
hole solutions are also obtained. 

PACS numbers: 04.20. J, 04.60.+n, 03.65. Ge 

1 Introduction. The model 

Multidimensional cosmology (see, for example [1-21] and references therein) is a very inter- 
esting object of investigations both from physical and mathematical points of view. 

Last decade the interest in multidimensional cosmology was stimulated mainly by Kaluza- 
Klein and superstring paradigmas [22,23]. The "realistic" multidimensional cosmological 
models appeared mainly in a context of some unifications theories. Certainly, it is quite 
natural to believe that the Entire Universe is multidimensional and we live in a some sort 
of a (3+1)- dimensional layer, that is Our Universe. Of course, at first stage we should try 
to understand the structure of our 3-dimensional crude (dense) matter and the formation of 
Our Universe. But it seems to be very likely that at some stage of our development it will 
be just impossible to describe our (3+l)-dimensional layer (Our Universe) out of touch with 
other (multidimensional) layers and domains. 

A large variety of multidimensional cosmological models is described by pseudo-Euclidean 
Toda-like systems [19] (see formula (1.10) below). These systems are not well studied yet. 
We note, that the Euclidean Toda-like systems are more or less well studied [24-28] (at least 
for certain sets of parameters, associated with finite-dimensional Lie algebras or affine Lie 
algebras). There is also a criterion of integrability by quadrature (algebraic integrability) 
for these (Euclidean) systems established by Adler and van Moerbeke [28]. Nevertheless, 
there are some indications that cosmological models may contain rather rich mathematical 
structures. For example, a self-dual reduction of the Bianchi-IX cosmology [29] lead us to 
the Halphen system of ordinary differential equations [30]. This system may be integrated in 
terms of modular forms [31] and is connected with a certain integrable reduction of the self- 
dual Yang-Mills equation [32] (with the infinite-dimensional group SDif f SU(2)). Another 
example is connected with the Kaluza-Klein dyon solution from [33]. The field equations 
for a spherically-symmetric Kaluza-Klein dyon in 5-dimensions were reduced in [33] to an 

1 This is an extended version of the talk presented by the authors at the 8th Russian Gravitational 
Conference (Pushchino, 25-28 May 1993) 
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open (Euclidean) Toda lattice with three points. Certainly, this problem may be formulated 
in terms of an appropriate cosmological model described by a pseudo-Euclidean Toda- like 
Lagrangian. So, we lead to an interesting nontrivial example of an integrable cosmological 
model. 

In this paper we consider a cosmological model describing the evolution of n Einstein 
spaces in the presence of m-component perfect-fluid matter. The metric of the model 

n 

g = -exp[2-f(t)]dt ^ dt + Y,exp[2x i (t)]g {i) , (1.1) 

i=i 

is defined on the manifold 

M = Rx Mi x . . . x M n , (1.2) 
where the manifold Mj with the metric is an Einstein space of dimension iVj, i.e. 

R mi nAg {i) ] = av/SU (i.3) 

i — 1, . . . , n; n > 2. The energy-momentum tensor is adopted in the following form 

m 

n = Y.TN {a \ (i.4) 

a=l 

(T^ a) ) = diag(-p< a \t),d a \t)SZ\. . . J:\t)5Z n ). (1.5) 
a = 1, . . . , m, with the conservation law constraints imposed: 

VmT; W =0 (1.6) 
cc = 1,...,to — 1. The Einstein equations 

<~^R = ^ (1.7) 

(k 2 is gravitational constant) imply y mT^ = and consequently SJmT^^ = 0. 

We suppose that for any cn-th component of matter the pressures in all spaces are pro- 
portional to the density 

p ( t\t) = (l-hf\x(t)))P (a) (t), (1-8) 

where 

t i a) (x) = ~^ a) (x), (1.9) 

i — 1, . . . , n, where $^(x) is a smooth function on R n , a — 1, ... ,m. 

In Sec. 2 the Einstein equations (1.7) for the model are reduced to the equations of 
motion for some Lagrange system with the energy constraint E — imposed. When m = 1 
and all spaces are Ricci-flat (A 1 = in (1.3), % = 1, . . . ,n) such reduction was performed 
previously in [9]. 
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In Sec. 3 we consider the Einstein equations, when all spaces are Ricci-flat and h\ = 
const, i — 1, . . . , n, a — 1, . . . , m. In this case we deal with pseudo-Euclidean Toda-like 
system with the Lagrangian 

1 m 

L A = -G ij x i x j - ]T n 2 A {a) exp^fV), (1.10) 

2 a =l 

where sign(Gij) = (—,+,...,+) [14,15], ulf 1 = Nih^ and = const i = l,...,n, 
a = 1, . . . ,m. The Einstein equations are integrated in the following cases: 1) m — 1; 2) 
n = 2, m > 2, ^ 0, u {a) - u (1) = a = 1, . . . , m, where u 2 = G ij u iUj = 0, u ^ 0; 3) 

u {a) = b (a) u ^ u 2 < Q ? A (a) > 0, « = 1 , . . . , m. 

In Sec. 4 the Wheeler-DeWitt (WDW) equation for the model is considered. When all 
spaces are Ricci-flat the WDW equation is integrated in all listed above cases. For m = 1 the 
solution of the WDW equation, satisfying so-called quantum wormhole boundary conditions 
[34], are obtained. 



2 The equations of motion 

The non-zero components of the Ricci-tensor for the metric (1.1) are following 

n 

J Roo = -E^[^-7^ + (^) 2 ], (2-1) 

i=l 

n 

Rm tni = gVjX* + exp(2x l - 2 7 )(x* + x\Y. - 7))], (2.2) 

i=l 

i — 1, . . . , n. 
We put 

n 

7 = 7o = £A^ (2.3) 

i=i 

in (1.1) (the harmonic time is used). Then it follows from (2.1) and (2.2) that the Einstein 
equations (1.7) for the metric (1.1) with 7 from (2.3) and the energy-momentum tensor from 
(1.4), (1.5) are equivalent to the following set of equations 

1 TO 

-G^iP + V c + k 2 J2 P (a) exp(2 7o ) = 0, (2.4) 

1 a=l 

m n 

X 1 + x l exp(2x i - 2 7o ) = k 2 exp(2^) £ [p^ + (D - 2)-\p^ - ]T N jP f\ (2.5) 

Q=l j = l 

i = 1, . . . , n. Here 

Gij = NAj - N l N J (2.6) 
are the components of the minisuperspace metric, 

1 n 

K = --E \ l N % cxp{-2x* + 2 7o ) (2.7) 
1 i=i 
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is the potential and D = dimM = 1 + E™=i 

The conservation law constraint (1.6) for a G {1, ...,m} reads 

n 

+ J- NitfipW + pi a) ) = 0. (2.8) 
i=i 

We impose the conditions of state in the form (1.8), (1.9). Then eq. (2.8) gives 

pi°)(t) = A^exp[-2N lX l (t) + & a \x(t))], (2.9) 
where = const and eqs. (2.4), (2.5) may be written in the following manner 

1 m 

-Gijtfx* + V + k 2 J2 aH exp $ (Q) = 0, (2.10) 

2 a=l 

m 

X 1 + x l exp(2x i - 2 7o ) = -h 2 u\ a) A {a) exp(2x* - 2 7o + $ (a) ), (2.11) 

a=l 

2 = 1, . . . , n. In (2.11) we denote 

= NM a) = dM a \ u\ a) = G*uf\ (2.12) 

where [15] 



b 13 1 

are the components of the matrix inverse to the matrix (Gy) (2.6). 

It is not difficult to verify that equations (2.11) are equivalent to the Lagrange equations 
for the Lagrangian 

L = l-Gijx*^ - V (2.14) 

where 

m 

V = V(x) = V c (x) + K 2 A^esxp[^ a) (x)]. (2.15) 

a=l 

Eq. (2.10) is the zero-energy constraint 

E = ]-G ij x i x j + V = 0. (2.16) 

Remark 1. In terms of 1-forms = u^dx 1 , the relations (1.9) read: = d& a \ 
a — 1, ... ,m. In this case 

du (a) =0, (2.17) 

a = l,...,m. The set of eqs. (2.17) (on R n ) is equivalent to (1.9). An open problem 
is to generalize the considered here formalism for the following cases: a) du^ ^ for 
some a G {1, . . . , m}; b) du^ = for all a — 1, . . . , m, but are defined on an open 
submanifold Q G R n with the non-trivial cohomology group iJ x (f2, R) ^ 0. 
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Using eqs. (2.1) and (2.2), it is not difficult to verify that the Einstein equations (1.7) 
for the metric (1.1) and the energy- momentum tensor from (1.4), (1.5), (1.8), (1.9) are 
equivalent to the Lagrange equations for the following degenerate Lagrangian (see also [15]) 

L = - exp(-7 + -f (x))G ij x i x j - exp(7 - 7o(a;))V(a;) (2.18) 
(L = L(7, x, x)). Fixing the gauge 

7 = 7o(a0-2/(a;), (2.19) 
where f — f(x) is a smooth function on K n , we get the Lagrangian 

L f = ± vx ] H2fi.r))(; il r'.h' - exp(-2f(x))V(x). (2.20) 

For / = we have the harmonic-time gauge (2.3). The set of Lagrange equations for the 
Lagrangian (2.18) (or equivalently the set of the Einstein equations) with 7 from (2.19) 
is equivalent to the set of Lagrange equations for the Lagrangian (2.20) with the energy 
constraint imposed 

E f = l - exp(2/(z))G^V + exp(-2f(x))V(x) = 0. (2.21) 

Remark 2. We remind that the action of the relativistic particle of mass m, moving in 
the pseudo- Euclidean background space with the metric Gij(x) has the following form 

dr[G y (x(r))^ - — e(r)], (2.22) 

where e = e(r) is 1-bein. Comparing (2.18) and (2.22), we find that for V(x) > the cosmo- 
logical model (2.18) is equivalent to the model of relativistic particle with the mass m — 1, 
moving in the conformally-flat (pseudo-Euclidean) space with the metric Gij(x) = 2V(x)Gij 
(see also [37]). In this case e = 2V(x) expi^j — jo(x)) . For V(x) < we have a tachyon. 
The problem may be also reformulated in terms geodesic-flow problem for conformally-flat 
metric (this follows from (2.22) or from a more general scheme [43]). 

3 Classical solutions 

Now, we consider the following case: A* = (all spaces are Ricci-flat), uf^ = Nih^ = const, 
i = 1, ... Then V c = and we put $^ = uf^x 1 in (2.15). In this case the Lagrangian 
(2.14) has the form (1.10). 

Remark 3. The curvature induced term V c (2.7) may be generated in the framework of 
the model with the Ricci-flat spaces Mj by the addition of n new components of the perfect 
fluid with uf ] = 2Ni - 2b\ and k 2 A^ = -\ k N k /2, i,k = 1, ... ,71. The introduction of the 
cosmological constant A into the model is equaivalent to the addition of a new component 
with uf +1) = 2Ni and k 2 A^ = A. 
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3.1 One-component matter 

We consider the case m—1, A^> = A ^ 0. We denote hf~p = hi, u~p —m — iVj/ij. 
We remind [14, 15] that the minisuperspace metric 

G = G ij dx i ®dx i (3.1) 

has pseudo-Euclidean signature (— , +,..., +), i.e. there exist a linear transformation 

z a = V?x\ (3.2) 

diagonalizing the minisuperpace metric (3.1) 

n-l 

G = r] ab dz a <g> dz b = -dz° ®dz° + Y, dz i <g> dz\ (3.3) 

i=i 

where 

(% 6 ) = fa 06 ) = dia ff (-l, +1, . . . , +1), (3.4) 

a, b — 0, . . . , n — 1. 

Proposition 1. For any -u = (tij) G i? n , a^O, there exists a (nondegenegate) nxn matrix 
(Vf) such that 

^V^/ = Gij (3.5) 

and a) V t ° = Ui/V^t 2 , for w 2 < 0; b) V> = m/Vv?, for u 2 > 0; c) V? + V* = u u for w 2 = 0; 
Here and below (u = (ui) = (iVj/ij)) 

n in 

M 2 = = G ij UiUj = Y / N t (h i ) 2 + — — QTA^) 2 . (3.6) 

i=l i=l 

(We note that in notations of [14] w 2 = A(/i)/(2 - £>).) 

This proposition follows from the fact that < u,v >= G^UiVj is bilinear symmetric 
2- form of signature (—,+,...,+) and the following quite obvious 

Proposition 2. Let v & E = R n , n > 2, and <.,.>: E x E — > i? is a bilinear symmetric 
2-form of signature (— , +,..., +). Then there exists a basis v°, . . . ) v n ~ 1 in E, such that 
< v a , f b >= r] afc and a) v — v°, b) v — v 1 , c) v = v° + v 1 , in the cases: a) v 2 =< i>, t> >= — 1, 
b) v 2 = 1, c) f 2 = respectively. 

Let -u 7^ 0. In z = (^-coordinates (3.2) with the matrix (Vf) from the Proposition 1 
the Lagrangian (2.14) has the following form 

1 1 n_1 1 

L A = - Vab z a z b -V A = --(i ) 2 + £ - V A , (3.7) 



where 



V A = K 2 Aexp(2qz°), u 2 < 0, (3.8) 
= /t 2 Aexp(2gz 1 ), m 2 >0, (3.9) 
= n 2 Aexp(z° + z 1 ), u 2 = 0, (3.10) 
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is the potential (2.15). Here we denote 

2q=^P\. (3.11) 
The Lagrange equations for the Lagrangian (3.7) 

z a = -r] ab d b V A (3.12) 

with the energy constraint (2.16) 

E A = \v ab z a z h + V A = 0, (3.13) 

can be easily solved. We present the solutions, 
a) For u 2 < 

z^pH + q 1 , i = l,...,ra-l, (3.14) 

2qz° = y(t), (3.15) 

where p % and q % are constants and 

y{t)= \n[C/Dsmh 2 dVC(t-t ))],C ^0,D>0, (3.16) 

ln[4/£>(f-i ) 2 ], C = 0,D>0, (3.17) 

= ]n[-C/Dcosh 2 (-VC(t-t ))],C>0,D<0, (3.18) 

Here t Q is an arbitrary constant, D = —2u 2 k 2 A, C = —u 2 (p) 2 and (p) 2 = J27=i(P 1 ) 2 - 



b) For u > we have 



z^pH + q*, i = 0,2,...,n-l, (3.19) 
2qz 1 = y(t), (3.20) 



with (p) 2 = (p ) 2 - E^V) 2 in (3.15)-(3.18). 
c) m 2 = 0, ii ^ 0. In this case 

z^pH + q 1 , i = 2,...,n-l, (3.21) 

2 + = 2° + = p + t + g + , (3.22) 

z~ = z° - z 1 =p - i + 9~ + k 2 Az(£), (3.23) 

where for p + ^ 

=2(p + )- 2 exp(p + t + g+), P + P~ = (p) 2 (3.24) 
(p~ = for n = 2) and for p + = 

z(t) = t 2 expg+ {p) 2 + 2K, 2 Aexpq + = 0. (3.25) 

Here (p) 2 = E^V) 2 - 
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For u = we have 

z a =p a t + q a , a = 0,...,ra-l, (3.26) 
|wV + k 2 A = 0. (3.27) 

Kasner-like parametrization. Here we consider the case it 2 < 0, 4 / 0. For C = 

—u 2 {p) 2 > we reparametrize the time variable 

r ln exp(Vg( t - f0 )) + yi 

V? exp(VC(t - t„)) - 

where 

e = A/|A| = ±l, T= (2//t 2 |A|| M 2 |) 1/2 . (3.29) 
We introduce new (Kasner-like) parameters 



a* = -Wlf I ^ -u\$)\ (3.30) 

where (V a J ) = (V/ 1 )" 1 and the summation parameter s runs: s = 1, . . . ,n — 1. Then, due to 
relations (3.2), (3.5), (3.14)-(3.16), (3.18) and Proposition 1 we get the following expression 
for the metric (1.1) [40] 

9 = ~{f{{ai{r)f Ni - Ui )dT ®dr + ± a 2 (r)g^, (3.31) 
i=i i=i 

where 

a t (r) = A[ ^h(r^/T) ]2MVu2[ tanh(r^/2T) r ^ ^ 

i = 1, . . . , n; Ai > are constants and the parameters a 1 satisfy the relations 

Uit** = 0, (3.33) 
GijaW = -4/u 2 (3.34) 

(see Proposition 1 and (3.30)). For the density (2.15) we have 

p(r) = Af[(a l (r)r- 2 ^. (3.35) 
i=i 

We note, that (p) 2 = 2k 2 \A\ U?=i A? 

For A > we have an exceptional solution (3.31), (3.33), (3.34) with the scale factors 

ai{r) = A i exp(±2u i r/u 2 T), (3.36) 

A\ > 0, i — 1, . . . , n. This solution correspond to C = case (3.17). 

Remark 4. In [19] the Einstein equations (2.10), (2.11) were solved for = 0, 

a = 1, . . . , m, A 1 7^ 0, A* = 0, % > 1. The solutions [19] may be also obtained from the 
formulas (3.31)-(3.34). We note that the spherically-symmetric analogue of the solution 
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[19] was considered in [36] (the case d — 2 was considered previously in [35]). There exists 
an interesting special case of the solutions [35, 36]. It is the n-time generalization of the 
Schwarzschild solution 

g = -[(1- ±) A ] ab dt a ® dt b + (1 - ±)- spA dR®dR+{l- ^f'^lPdVL 2 , where L ^ and 
A = (A a b) is symmetric n x n matrix, satisfying the relation sp(A 2 ) + (spA) 2 = 2. 

We consider this solution in a separate publication. 

3.2 Two spaces with m-component matter 

Here we consider the following case: n = 2, m > 2, A {a) ^ 0, 



(«)_ u (i) = &(«) u (3.37) 



u 

a = 1, . . . , m, where u 2 = 0, u ^ and are constants. 

In ^-coordinates (3.2), where the matrix (V?) satisfies the Proposition 1 (see the case c) 
u 2 = 0) we have 

z + = z ° + z 1 = (y° + Vl)^ = Ui x\ (3.38) 
$M = Mj ( V = + a-z~, (3.39) 

where 2a + = — < >, 2a_ = — < v},u >, and u* = (u*) is defined by the relation : 

u*x l = z~ (or equivalently < u*,u* >= 0, < u*,u >= —2). 
Due to (3.37)-(3.39) the potential in (1.10) is factorized 

V = V + (z + )V-(z-), (3.40) 

where 

m 

V+(z + ) = exp(a + z + )(K 2 A^ + £ k 2 A^ exp(b {a) z + ), (3.41) 

k=2 

V-(z~) =exp(a_z-). (3.42) 
Let > 0, a = 1, . . . , m, . We consider the /-gauge (2.19) with 

F = e 2f = V. (3.43) 
In this gauge the Lagrangian (2.20) reads 

Lf = -\v + {z + )z + V_{z-)z- - 1. (3.44) 



In the variables 



w ± 



w ± {z ± ) = f dxV±(x) (3.45) 

Jzo 

the Lagrangian (3.44) has rather simple form 

L f = -\w + w- - 1. (3.46) 
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The equations of motion for (3.46) give 

w ± (t) =p ± t + q ± . (3.47) 
The parameters p ± satisfy the energy constraint 

2E f = -p + p~ + 2 = 0. (3.48) 

Remark 5. It is interesting to note that the so-called D-dimensional Schwarzschild- 
deSitter solution [44,45] may be obtained from the considered here cosmological solution 
with n = m = 2 and JV X = 1, N 2 = D - 2. 

3.3. n spaces with m component matter 

Now we consider the simplest case of the multicomponent matter. We put in (1.10) 
n > 2, > 0, = b^u, u 2 < 0, where are constants, a — 1, . . . , m. 

In ^-coordinates (3.2), corresponding to the case a) from the Proposition 1, the La- 
grangian (1.10) has the form (3.7) with the potential 

m 

V A = V A (z°) = ]T K 2 A ia) exp(2qb {a) z°), (3.49) 
i=i 

where q is defined in (3.11) (A = (A^)). The solutions of the equations (3.12) and (3.13) 
are expressed by the formula (3.14) and the following relation 

/ dx[2S + 2V A (x)}- 1/2 = ±(t-t ), (3.50) 

where 28 = E^iiP*) 2 , and c ,t are constants. 



4 Quantum solutions 

The WDW equation for the model in harmonic time gauge (2.3) reads as follows: 

(—^■+^ = 0, (4.1) 

where \I/ = ty(x) is "the wave function of the Universe", V is the potential (2.15), <9j = d/dx 1 
and G tj are defined in (2.13). The relation (4.1) is a result of a trivial quantization of the 
zero energy constraint (2.16), written in the form \iE = 0. Here fj, is a fundamental quantum 
parameter of the theory. 

In /-gauge (2.19) the WDW equation should be written in the conformally covariant 
form [15, 37] (such form of the WDW equation was discussed earlier by Misner [38]) 

(-^-A[e 2f G] + —R[e 2f G] + e~ 2f fiV)V f = 0, (4.2) 

where A[G] and R[G] are the Laplace-Beltrami operator and the scalar curvature of G 
respectively, a n = (n — 2)/8(n — 1) and 

* 7 = exp[(2 - n)//2]tf. (4.3) 
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Without loss of generality we put pi — 1 below. 
4.1 One-component matter 

Here we find the quantum analogues of the classical solutions from 3.1, i.e. we integrate 
the WDW equation 

<4"°^+ iy *= ' < 4 - 4 » 

with the potential (3.8)-(3.10). We note, that the WDW equation for 1-component model 
with n Ricci-flat spaces was considered previously in [17]. 

a) u 2 < 0. In this case the WDW equation (4.4) reads 

te-Sss 1 ^ 1 *' 1 (4 ' 5) 

We are seeking solutions of (4.5) in the following form 

V(z) = exp(ipz)<Z>(z°), (4.6) 

where p = (pi, . . . ,p n -i) is a constant vector (generally from C n ~ v ), z — (z 1 , . . . , z n_1 ) , pz = 
YaZ\Vi z1 1 ■ The substitution of (4.6) into (4.5) gives 

[-(A) 2 - 2^exp(2^ )]$ = 2£$, (4.7) 

where 2S = Y%=iPi- Solving (4.7), we get two linearly independent solutions 

$( z °) = B v (y/-2K?Aq- 1 e q '°), (4.8) 
where v = iV2£/q = i\p\/q, and B u = I u , K v is modified Bessel function. We note, that 

1 n 

v = exp qz° = exp(-UiX l ) = J| a" i/2 (4.9) 

is a natural scale factor for the model (aj = e x '). 

The general solution of eq. (4.5) has the following form 

*(*)= E / ' d^pCsimfiz), (4.10) 

B=I,K J 

where 

*j!(z) = e^B m/q (V^Aq^ z °), (4.11) 

and functions Cb (B = I,K) belong to an appropriate class of (generalized) functions. 

b) u 2 > 0. In this case the WDW equation (4.4) reads 

8 8 8 8 ™~ 1 8 8 

^TO + W ~ g a? a? + ^«pM1* = o. (4,2) 
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An analogous consideration in this case gives the general solution (4.10) with 

tff(z) = e^B^iV^Aq-^ ). (4.13) 



Reiep = (p ,p 2 , • • ■ ,Pn-i), 2 = (^ W . . , £ n x ), v(p) = iV2£/q, and 2£ = p 2 - Yl?J pi 
c) u 2 = for -u 7^ the WDW equation reads 

n-l a 

[-40 + 0_ + + J2(jri) 2 - 2« 2 Aexp(z + )]tf = 0, (4.14) 
i=i ° z 

where z ± = z° ± z 1 , d± = d/dz ± . The substitution 

^f(z) = exp(ipz)$(z + , z~), (4.15) 
with p = (p 2 , ■ ■ .,p n -i), z= (z 2 , . . . , z n ~ l ) entails 

[4<9+<9_ + 28 + 2K 2 Aexp(z + )]<5> = 0, (4.16) 
where 2£ = Yh=2 Pi ■ Introducing new variables u°, u 1 , where u° ± u 1 = u 41 and 

u + = 2Sz + + 2n 2 Aexp(z + ), vT = z~ (4.17) 
we get the Klein-Gordon equation for $ with m 2 = 1 

0"- ( s? ),+1)# - a (418) 

It is quite obvious how to write the general solution of (4.22). 

Quantum wormholes. In the case a) u 2 < for A < there exist so-called quantum 
wormhole solutions of the WDW equation [34]. We present here a continuous spectrum 
family of these solutions. The wave functions are following 

if X:H (z) = exp[-q'W-2K 2 Ae qz ° cosh(A - qzn)}. (4.19) 

where A G R and n is unit vector: (n) 2 = 1 {n G S n ~ r ). These solutions are related with the 
solutions (4.11) (with B = K) by the formula 

= - / dk^ qkn {z)e- ik \ (4.20) 

71 J-oo 

(such trick was suggested in [39], see also [20,41]). The solutions (4.19) satisfy the quantum 
wormhole boundary conditions (in terms of parameter v (4.9): i) the wave function is expo- 
nentially damped for large space geometries (v — > +oo); ii)the wave function is regular when 
the spatial geometry degenerates (v — > 0). 
We also note that the the functions 

( r 0\2 , / 1\2 

V m , n = HUx^HUx 1 ) exp[- 1 ) ^ [ ) ] (4.21) 
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where H m are Hermite polynomials, m — 0, 1, . . . „ 

x i = {2/qfl 2 {-2K 2 A) 1 l A exp{qz°/2)f{±qzn), ' i = 0, 1, 

(/°) f 1 ) = {sink, cosh) are also solutions of the WDW equation with the quantum worm- 
hole boundary conditions. (They are called discrete spectrum quantum wormholes.) We 
note that the special cases of the solutions (4.19), (4.21) for Ui = 2iVj (A-term case) and 
Ui = 2Ni — 28} (1-curvature case, A 1 ^ 0) were considered in [41] and [20] respectively. 

We also note that for b) u 2 > and A > there also exist quantum wormhole solutions. 
(In this case z° should be replaced by z 1 in (4.19), z is defined in 3.1 b) and n belongs to 
hypersphere.) 

4.2. Two spaces with m-component matter 

For the model from the subsection 3.2 the WDW equation (4.2) in the /-gauge (3.43) 
has the following form (/i — 1) 

( v +( J)fc+vi(JW +l) * =0 - (422) 

Indeed, for n = 2 we have A[e 2 ^G] = e~ 2 ^A[G] , <i2 = and tyf = (see (4.3). In w-variable 
w = (w ,^ 1 ), where w°±w 1 = w ± , where w^ 1 are defined in (3.45), we get the Klein-Gordon 
equation with m 2 = 2 

[( ^-(^ + 2]* = 0. (4.23) 
4.3 n-spaces with m component matter 

Here we present the solutions of the WDW equation (4.4) with the potential (3.49) , i.e. 
quantum analogues of the classical solutions from 3.3 are considered. 

Repeating all arguments from 4.1 (case a)), we get the general solution of (4.4) 

= W <f-Vc t ( P >;(4 (4.24) 

*=± J 

where 

Vt{z) =exp(ipS)^°), (4.25) 
* = ±, and $* = &fi(z°) are two linerly independent solutions of the equation 

[-(^) 2 -2^(z°)]$ = 2^<I>, (4.26) 

with the notations for S = £p, p, z from 4.1. a). 

We note, that for special values of parameters and in the potential (3.49) the 
equation (4.26) describes the quantum spin systems [42]. 



5 Concluding remarks 

In this paper we investigated the multidimensional cosmological model with n (n > 1) Ricci- 
flat spaces, filled by m-component perfect fluid. In some sense, this model may be considered 
as "universal" cosmological model: a lot of cosmological models may be obtained from it 
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under a suitable choice of parameters. This fact may be used for "Toda-like" classification of 
known exact cosmological (and spherically-symmetric) solutions of the Einstein equations. 
(We note, that the Bianchi-IX cosmological model is described by the "Toda-like" Lagrangian 
(1.10) with n = 3 and m = 6.) 

Here we integrated the Einstein and Wheeler-DeWitt equations for some sets of param- 
eters. But an open problem is the problem of integrability of the considered here model (at 
classical and quantum levels) for arbitrary values of the parameters m, n, Ni and u^. We 
hope to continue the investigation of this problem in forthcoming publications. 
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